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SECTION  1 


INTRODUCTION  AND  OVERVIEW 

1.1  BACKGROUND  AND  BASIS  OF  PRESENT  STUDY 

The  present  study  was  prompted  as  a result  of  various  com- 
munications and  publications  of  Cooper  (Reference  1 for  example) . 
He  correlated  the  close-in  horizon cal  ground  displacements  of 
cratering  experiments  to  the  apparent  crater  volumes  and  suggested 
that  his  results  were  a manifestation  of  a late  time,  essentially 
incompressible  flow  phenomena  that  should  be  explained. 

John  Lewis*  noted  the  existence  of  a large  body  of  printout 
data  corresponding  to  computer  simulations  of  the  early  cratering 
flow  that  occurs  before  material  strength  effects  appear.  He 
suggested  the  possibility  of  discovering  properties  of  these 
flows  which  could  provide  a basis  for  understanding  Cooper's 
observations  and  the  mechanics  of  cratering  and  ejecta  to  pro- 
vide scaling  rules  and  superior  prediction  methods.  This  is  the 
basis  of  the  computational,  experimental,  and  correlative  study 
that  is  reported  here. 

This  study  was  initiated  by  DNA  Contract  DNA001-72-C-0063 
and  a draft  final  report  was  submitted  in  February,  1973.  The 
results  at  that  time  were  premature  and  a DNA  document  number 


¥ 


Private  communication  from  Mr.  Lewis  of  DNA, 


presently  at 


RDA. 
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was  not  assigned.  The  efforts  were  continued  without  DNA  support 
and  a first-order  cratering  flow  model  was  developed  by  late  1973. 
This  model,  based  on  five  surface  burst  computations,  was  in 
satisfactory  agreement  with  limited  experimental  and  computed 
data.  DNA  support  returned  in  late  1973. 

The  primary  objective  of  the  present  study  was  to  assess  and 
confirm  the  first-order  flow  model.  This  was  accomplished  satis- 
factorily. The  second  objective  was  to  improve  and  extend  the 
model  to  provide  improved  prediction  methods  and  scaling  rules 
for  cratering,  displacements,  and  ejecta.  Substantial  progress 
was  made  toward  this  objective. 

Section  2 reports  the  steady- state  flows  that  were  observed 
in  10  surface  burst  cratering  calculations  and  develops  and  as- 
sesses the  first-order  flow  model.  A few  confirming  examples  are 
given.  Many  others  were  redundant  and  were  omitted.  No  serious 
discrepancies  were  found. 

Section  3 develops  and  applies  an  approximate  energy  balance 
model  to  freeze  the  cratering  flow  into  the  final  crater  size  and 
shape  by  employing  an  elastic-perfectly  plastic  material  assump- 
tion. Plasticene  clay  cratering  data  was  predicted  satisfactorily 
by  this  method.  Other  real  materials  appear  to  have  very  little 
shear  strength  along  the  flow  paths  that  lead  to  the  final  crater 
radii,  and  the  effect  of  gravity  is  significant  in  these  cases. 

Section  4 is  a potpourri  of  data  ind  considerations  that 
relate  to  advanced  models.  Because  the  final  model  does  not 
exist,  it  was  not  possible  to  organize  this  section  in  a satis- 
factory way.  A large  amount  of  data  of  unknown  significance  has 
been  omitted. 
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Section  5 presents  a summary  and  conclusion. 


Appendix  A reports  the  experimental  program  that  was  conduc- 
ted in  support  of  the  theoretical  and  correlative  efforts. 

Appendix  B gives  some  of  the  mathematical  details  that 
relate  to  the  energy  balance  method  of  Section  3. 

The  remainder  of  this  section  will  present  an  overview  of 
Sections  2 and  3. 

1.2  OBSERVATIONS  AND  MODELING  OF  STEADY-STATE  CRATERING  FLOW 

Surprisingly  simple  cratering  flows  were  observed  in  all 
cases  that  were  examined.  Steady-state  flows  in  the  earth 
(independent  of  time  at  a fixed  point  in  space)  were  set  up  in 
the  cratering  region  immediately  after  shock  wave  passage  when 
the  narrow  stress  pulse  had  relaxed  and  the  material  density 
approached  its  final  constant  value. 

Mass  continuity  expressed  as  V • pu  = 0 zould  then  be 
approximated  by  V • u = 0 corresponding  to  incompressible  flow. 
Then  a steady-state  approximation  of  one  velocity  component 
fixes  the  other,  A first  order  flow  model  was  constructed  using 
this  method. 

The  first  order  model  employs  two  constant  parameters?  an 
intensity  parameter,  a,  and  a shape  parameter,  Z,  to  approximate 
the  flow  everywhere  in  space  to  predict  average  cratering 
response. 
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1.3  TERMINATION  OF  THE  STEADY-STATE  FLOW 


The  steady-state  flow,  established  at  a time,  tQ,  shortly 
after  shock  passage,  was  observed  to  persist  for  a long  time 
which  could  be  orders  of  magnitude  larger  than  tQ  (in  the  case 
of  weak  soil  or  sand) . When  material  strength  or  gravity  effects 
become  significant,  the  flow  slows  markedly  and  stops  without 
appreciably  changing  direction  (Z  remains  constant) . The  final 
crater  depth  occurs  first.  The  process  of  slowing  and  stopping 
occurs  later  at  higher  angles  with  the  final  crater  radius 
occurring  last. 

The  method  that  emerged  to  predict  the  termination  of  the 
flow  employs  stream  tubes  defined  by  adjacent  steady-state  stream 
lines.  It  is  assumed  that  energy  balance  can  be  applied  to  the 
mass  in  a stream  tube  independent  of  the  others.  The  kinetic 
and  gravitational  energies  depend  on  the  flow  parameters  a and  Z. 
The  distortional  energy  is  more  complicated  and  requires  a ma- 
terial model.  The  distortional  energies  were  evaluated  by  as- 
suming a constant  shear  modulus,  u,  and  a constant  von  Mises 
yield  stress,  Y.  The  resulting  expressions  for  kinetic,  gravita- 
tional and  distortional  energies  were  then  used  to  approximate 
energy  conservation  in  the  stream  tube  with  the  flow  stopping 
when  the  kinetic  energy  was  consumed. 

Plasticene  (children's  modeling  clay)  is  a reasonable  ap- 
proximation of  a constant  yield  stress  material  (except  for 
creep  and  temperature  dependency)  and  it  was  employed  in  an 
experimental  cratering  program  with  three  different  explosive 
heights  of  burst,  HOB.  The  corresponding  values  of  u,  z 
obtained  by  short  computer  si/  ulations  resulted  in  good  agree- 
ment between  experimental  and  first-order  model  crater  dimensions. 
However,  few  real  earth  media  can  be  approximated  properly  with 
an  elastic-perfectly  plastic  model  and  a more  general  method  is 
required. 
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1.4  MATERIAL  STRENGTH  PROBLEM  AND  SOLUTION 


Though  it  is  obvious  that  most  earth  materials  cannot  be 
molded  satisfactorily  by  a constant  yield  stress  assumption,  it 
is  not  known  what  models  are  appropriate  for  the  shock  condi- 
tioned, severely  strained  and  failing  material  in  the  cratering 
flow.  Stress-strain  paths  encountered  in  laboratory  triax  tests 
are  far  from  those  desired.  It  is  thought  that  the  state  of  the 
art  for  strength  modeling  of  earth  media  is  not  satisfactory  and 
a new  approach  must  be  taken  to  establish  effective  material 
strengths  for  the  cratering  process. 

A straightforward  and  attractive  solution  is  to  invert  the 
prediction  formulas  and  use  experimental  crater  depths  \nd  radii, 

D and  RAD,  to  establish  different  effective  yield  stresses,  YQ 
and  to  account  for  the  different  stress-strain  paths  en- 

countered in  the  crater  depth  and  crater  radius  flows.  A table 
of  values  of  YD  and  Yp^  for  a spectrum  of  materials  can  be  estab- 
lished in  this  way  to  provide  a means  of  estimating  the  effective 
strengths  of  similar  materials.  The  predictive  methods  then 
correspond  to  scaling  rules  where  only  the  u,Z  values  of  the  new 
case  are  required. 
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SECTION  2 

STEADY-STATE  FLOWS,  FIRST-ORDER  MODELS 


2.1  DISCUSSION 

The  computer  simulations  of  Table  1 provided  the  bulk  of  the 
data  that  led  to  the  flow  models.  The  first  10  cases  are  surface 
bursts  computed  by  two-dimensional,  finite-difference  codes.  The 
eleventh  case  is  a one-dimensional,  spherical  computation  of 
PILEDRIVER,  a deeply  buried  nuclear  event.  It  is  included  to 
illustrate  a feature. 

Figure  1 shows  schematic  examples  of  a typical  two- 
dimensional  case  showing  an  original  grid  configuration  in  the 
crater  region  and  its  distortion  and  velocity  field  at  an 
intermediate  time  (after  material  densities  return  to  a near 
constant  value  and  before  material  strength  effects  appear) . 

Referring  to  the  spherica]  coordinate  system  of  Figure  lc, 
it  was  observed  that  the  radial  velocity  was  approximately 
independent  of  angle  in  the  region  below  the  original  ground 
surface  (9  < tf/2)  . It  was  also  seen  that  this  velocity  could 
be  approximated  by 


R = a/RZ  (1) 

Applying  the  incompressibility  approximation,  ,V  • u = 0,  the 
tangential  velocity  component  is 

R9  = R(Z  - 2)  sin0/(l  + cos6)  (2) 
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TABLE  1 

COMPUTER  SIMULATIONS 


Case 

Name 

Simulation 

1 

ELK  64 

JOHNIE  BOY,  500- ton  nuclear  event.  DOB  = 1. 
feet  in  NTS  Area  18  alluvium 

2 

ELK  67 

Z-C6,  Plasticene  crater  experiment  scaled  to 
1 cm  radius  HE  sphere.  HOB/R™  = 0 (half- 
buried  sphere) 

3 

ELK  66 

Z-C2,  same  as  above  except  HOB/R™  = 1 
(tangent  sphere) 

4 

ELK  69 

Z-C7,  same  as  above  except  HOB/R^  = 2 

5 

ELK  65 

CACTUS,  17-kt  nuclear  surface  burst  over 
5.5  feet  of  coral  sand  overlaying  beachrock 

6 

NOT  DISCUSSED 

7 ELK  73  5-Mt  nuclear  surface  burst  over  layered  site 

of  shale  and  limestone 

8 ELK  76  Same  as  above  with  slightly  enhanced  energy 

due  to  neutron  deposition 

9 ELK  70 (L)  Surface  pressure  loading  function  P(R,t)  over 

coral  sand  layers  with  low  porosity. 

10  ELK  75(H)  Same  as  above  with  high  porosity  layers 

11  PILED RIVER  PILEDRIVER,  61-kt  nuclear  source  buried  in 

granite  at  121-bar  overburden  depth. 
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The  particle  motion  in  this  velocity  field  is 


RZ+1(t)  - RZ+1 


RZ+1(t> 


= (Z+l)a  (t  - 


to> 


(3) 


l - cose  (t) 

1 - COS0 

o 


[R(t)/Ro]z_2 


(4) 


where  rQ,  dQ  are  the  particle  coordinates  at  tQ.  In  most  cases 
of  interest  the  onset  time  is  small  and  the  displacement  accumu- 
lated curing  shock  passage  is  negligible.  Then  tQ  % 0 and  Rq 
and  0Q  can  be  approximated  by  the  undisturbed  coordinates  of  the 
particle.  Rc(t)  is  defined  as  the  radius  of  a hypothetical 
particle  that  originated  at  Rq  = 0.  It  corresponds  to  the  radius 
of  a crater  wall  point  when  no  mass  has  been  lost  by  vaporization 
or  excavation  to  implant  the  explosive  source. 


The  steady  state  equations  (1-4)  were  compared  to  Table  1 
cases  and  good  agreement  was  observed.  A few  examples  will  be 
presented. 


2.1.1  Radial  Velocity.  Equation  (1)  was  a reasonable 
approximation  of  the  steady  state  radial  velocity  that  was 
observed  in  all  cases.  Two  cases  are  reported  as  examples. 

Case  11,  PILEDRIVER,  was  spherically  symmetric  correspond- 
ing to  Z = 2.  The  computed  radial  velocity  R versus  R at  various 
times  is  displayed  in  Figure  2.  The  narrow  propagating  peak 
corresponding  to  the  shock  wave  can  be  seen.  The  transition  to 
approximate  steady  state  flow  occurs  immediately  behind  the 
narrow  shock  with  Equation  (1)  being  an  excellent  model  until 
slightly  after  100  msec  when  material  strength  and  overburden 
effects  become  significant. 
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A one-dimensional  representation  of  R versus  R can  also  be 
constructed  for  the  two-dimensional  flows  in  the  manner  given 
next. 


Referring  to  Figure  1,  let  S..  be  any  of  the  original 
hemispherical  surfaces  that  surround  ground  zero.  The  volume 
below  the  original  ground  plane  and  inside  is 


2tt 

3 

tt/2 

r*  <e,t) 

C 

vi(t) 

/ 

1 R?(0,t)dRi  sinQ  d0 

(5) 

Jo 

o 

Similarly, 

the  rate  of 

change  of  volume  V. (t)  can  be 

evaluated . 

The  corresponding 

average  radius  and  radial  velocity 

are  then 

defined. 

/a  \1/3 

' (' T*) 

(6) 

• 

R. 

1 

= V±/  (2ttR2) 

(7) 

In  Case  1,  JOHNIE  BOY,  the  original  hemispherical  grid 
lines  with  Rq  (meters)  = 3.6,  3.8,  4.0,  4.2,  were  examined 

at  various  intermediate  times  after  shock  passage  and  values  of 
and  R^  were  constructed.  The  results  are  displayed  in  Figure  3. 
The  notations  3.6,  3.8,  etc,,  refer  to  the  value  of  Rq.  Points 
ahead  of  and  in  the  shock  wave  are  omitted  for  clarity.  The 
agreement  with  Equation  (1)  is  good  until  shortly  after  100  msec. 

This  method  of  constructing  average  values  of  a and  Z was 
also  used  for  Cases  2 to  5.  The  results  are  listed  in  Table  2. 
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TABLE  2 


AVERAGE  VALUES  OF  a 

AND 

Z CONSTRUCTED  BY 

VOLUME  METHOD 

Case 

Name 

, Z+l  . . 

a (cm  /\i  sec) 

Z 

1 

500-ton  JOHNIE 

BOY 

0.6422  x 106 

2.71 

2 

Plasticene  HOB 

= 

0 

0.563 

3.0 

3 

Plasticene  HOB 

= 

1 

0.151 

3.0 

4 

Plasticene  HOB 

= 

2 

0.0405 

2.7 

5 

17-kt  CACTUS 

2.1  x 106 

2.75 

11 

61-kt  PILEDRIVER 

3.8  x 105 

2.0 

2.1.2  Velocity  Field.  Space- averaged  values  of  a and  Z 
used  in  Equations  (1)  and  (2)  correspond  to  the  first-order 
model  of  the  steady-state  flew  field.  This  model  was  in  fair 
agreement  with  the  observed  data  in  all  cases  with  the  largest 
discrepancies  near  the  surface.  Two  cases  are  given  next  as 
examples. 

Figure  4 displays  Case  1,  JOHNIE  BOY  velocity  vectors  at 
selected  points  and  selected  times  as  indicated  by  the  symbols 
and  solid  arrows.  Only  a few  of  the  many  points  are  shown.  The 
dashed  arrows  are  the  model  velocities  constructed  from  Equa- 
tions (1)  and  (2)  and  the  average  values  of  a and  Z from  Table  2. 
Though  the  overall  agreement  is  satisfactory,  there  is  a definite 
tendency  to  the  discrepancies.  The  model  vectors  point  above 
the  computed  vectors  at  low  angles  and  below  the  computed  vectors 
at  high  angles.  Z,  the  flow  shape  parameter,  has  an  angular 
dependency.  Its  radial  and  time  dependency,  if  present,  is  not 
obvious. 


17 


% 

V 

°s* 

Computer  velocity 

Model  velocity 

\ 

% 

Symbol 

Time  Velocity  scale 

(msec)  (m/sec) 

Horizontal  range  (m) 


Figure  4 Case  1.  Equation  (2)  velocities  (dashed  arrows)  compared 

to  computed  velocities  (solid  arrows)  at  selected 
points  and  times. 


The  layered  Case  5,  CACTUS,  is  chosen  for  a similar  display. 
The  velocity  fields  are  shown  in  Figure  5.  The  discrepancies  are 
similar  to  those  of  the  unlayered  case  discussed  above  and  the 
same  remarks  apply. 

2.2  STREAM  LINES  AND  DISPLACEMENT  FIELD 

The  stream  lines  modeled  by  Equations  (3)  and  (4)  are  shown 
in  Figure  6 for  Z = 2.71.  This  model  is  clearly  non-physical  at 
high  angles  where  the  stream  lines  become  vertical  and  then  curve 
back  to  the  symmetry  axis. 

The  time  dependent  displacements  of  selected  points  at 
various  times  are  displayed  in  Figure  7 for  Table  2,  JOHNIE  BOY 
values  of  a and  Z.  The  very  thin  ejecta  lip  shape  cannot  be 
expected  to  occur  for  a variety  of  physical  reasons.  The  smooth- 
er portions  of  the  surface  upthrust  were  in  satisfactory  agree- 
ment with  the  computer- simulated  counterpart. 

2.2.1  Trajectories.  The  steady-state  velocity  vector 
fields  prescribe  the  particle  trajectories  and  a satisfactory 
model  of  one  necessarily  is  satisfactory  for  the  other.  One 
trajectory  example  will  be  displayed,  however,  to  validate  an 
approximation  used  later.  The  approximation  is  that  the  stream 
tubes,  defined  by  adjacent  stream  lines,  do  not  interact  so  that 
the  flow  stopping  in  one  does  not  change  the  direction  of  the 
others.  This  assumption  decouples  the  flow  that  defines  final 
crater  depth  (which  occurs  first)  from  that  of  the  final  crater 
radius  (which  can  occur  much  later) • 

The  time-dependent  crater  profile  of  Case  1,  JOHNIE  BOY,  is 
shown  at  various  times  in  Figure  8.  The  final  crater  depth  is 
reached  in  about  200  msec  with  the  final  radius  occurring  at 
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5.  Equation  (2)  velocities  (dashed  arrow)  compared 
to  computed  velocities  (solid  arrows)  at  selected 
points  at  5 msec. 


(c)  c • 200  msec 


(d)  t » 790  msec 


Figure  7 


Progressive  stages  of  cratering  by  Equation  (3) 
model  (see  text) . 


about  800  msec.  The  trajectories  of  points  on  the  crater  wall 
and  inside  the  crater  wall  were  examined  and  were  found  to  agree 
satisfactorily  to  the  Equations  (3)  and  (4)  model  when  the 
constant  a,  Z values  of  Table  2 were  employed.  Figure  9 displays 
this  comparison  for  particles  that  originated  at  a spherical 
radius  of  10  meters.  As  expected,  the  discrepancies  are  consis- 
tent with  those  of  the  velocity  fields  of  Figure  4.  The  computed 
particle  radii  at  200  msec  are  in  good  agreement  with  the  model. 
Material  strength  effects  then  appear  and  the  computed  motion 
slows  and  stops,  first  in  the  crater  depth  region  (0  - 0)  and 
later  at  higher  angles.  The  particles  do  not  depart  substantially 
from  their  steady-state  streamlines  during  the  termination  of  the 
flow. 


2 . 3 EJECTA 

The  steady-state  flow  approximation  applied  at  the  surface 
(0  = tt/2)  can  be  used  to  model  the  process  of  mass  ejection 
through  this  surface.  Figures  10a  and  10b  are  schematics  showing 
the  surface  velocity  vectors  and  the  vertical  components  at  a 
time  t^.  Figures  10c  and  lOd  are  corresponding  displays  at  a 
later  time,  t2.  The  flows  are  identical  in  regions  where  mass 
is  present. 

It  is  convenient  to  define  a vertical  velocity  cutoff, 
xcut'  to  identify  mass  that  is  ejected  with  vertical  velocities 
equal  to  or  greater  than  a prescribed  value.  The  following 
definitions  and  expressions  are  used: 
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(C)  Surface  velocities  at  t2 


(B)  Surface  vertical  velocity 
components  at  t-j 


(D)  Surface  vertical  velocity 
components  at  t2 


Figure  10  Schematics  of  surface  ejection  sequence. 


26 


I 

i 

E 


i 


p 

X 


y 

4> 

ccut 


cut 

Rc(t) 

RAD 


average  density  during  ejection 
2 

(Z  - 2)  q/R  = steady-state  vertical  velocity  at  6 = tt/2 
2 

o/R  = steady- state  horizontal  velocity  at  0 = tt/2 
tan  ^ (Z  - 2)  •=  ejection  angle  measured  from  horizontal 
reference  cutoff  for  vertical  velocity 

• 1 /7 

(Z  - 2) a/xcut  ' = radius  at  which  xcut  occurs 

[(Z  + l)atJ1/^^Z+1^  = crater  radius  at  time  t 
final  crater  radius,  known  or  estimated 
RADZ+1/[  (Z  + 1)  a]  =*  time  of  final  crater 


The  steady-state  velocity  expressions  are  not  valid  when  R 

c 

approaches  RAD  (or  t approaches  tf) . In  particular,  the 
velocities  are  0 at  t = tf  and  ejection  stops.  However,  the 
mass  ejected  with  velocities  x > xcut  is  unaffected  by  the  final 
crater  size  or  final  crater  time  provided  that  Rcut  < RAD.  This 
is  the  reason  why  Rcut  and  xcut  are  particularly  useful  parameters 


A variety  of  ejecta  expressions  can  be  derived  by  simple 
integrations.  Three  of  these  will  be  given  as  examples.  The 
explicit  time  dependencies  in  these  ejecta  expressions  are  re- 
placed by  the  more  convenient  parameter,  Rc(t) . Similarly,  xcut 
is  replaced  by  Rcu+.* 


M(Rcut'Rc> 
with  x > x 


= rate  of  mass  ejection  at  time 
cut 


t 


M(icut't) 


M(Rcut'Rc> 


X > X 


cut 


= mass  ejected  up  to  time  t with 


M(*cut'  V = M^Rcut,RAD)  = total  mass  ejected  with  x > xcut 


The  results  are 


M(xcut't) 


M(xcut'tf) 


Rc  < Rcut 


Z-2 


R < R . 
c cut 


cut 


< RAD 


(8) 


(9) 


(10) 


These  expressions  do  not  discriminate  between  the  surface  upthrust 
mass  (which  stays  in  contact  with  the  earth  at  all  times)  and  the 
true  ejecta  which  is  lofted  on  ballistic  trajectories.  However, 
it  can  be  shown  that  Equation  (10)  essentially  predicts  the  lofted 
mass  when  Rcut  = RAD.  The  upthrust  mass  is  then  the  difference 
between  the  total  mass  displaced  by  the  crater  and  the  lofted 
mass. 


Only  Case  1 (JOHNIE  BOY'  and  Case  8 (ELK  76)  were  computed 
to  final  crater  time  to  provide  enough  ejecta  data  for  sensitive 
comparison  to  Equation  (10) . The  agreements  were  fair  when 
volume  averaged  a and  Z values  were  used  (first-order  model)  and 
would  be  improved  if  surface  average  values  of  a and  Z were  used. 

The  comparisons  for  JOHNIE  BOY,  with  the  a and  Z values  of 
Table  2,  will  be  displayed.  Recalling  previous  remarks,  the 
JOHNIE  BOY  velocity  vectors  were  above  those  of  the  first-order 
model  at  high  angles  (see  Figure  4) . Then  at  a cutoff  range. 
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Rcut'  the  ratio  x/y  will  not  equal  the  volume  average  value 
Z - 2.  Therefore,  the  expressions  M(xcut*tf)  and  M( (Z-2)ycufc,tf) 
are  not  identical  and  it  is  instructive  to  compare  both  to  the 
computed  data. 


The  model  results  are  shown  as  the  dashed  line  in  Figure  11. 
The  circle  points  were  obtained  summing  the  computed  ejecta  mass 


with  the  constraint  x > x 


cut’ 


The  diamond  points  correspond  to 


the  sum  with  y(Z-2)  > xcut’  The  lar9e  discrepancies  that  occur 
above  velocities  of  about  10~^  may  be  due  to  spall  and  the  asso- 
ciated error  in  the  constant  density  assumption.  Otherwise  the  dis 
crepancies  between  the  model  and  either  of  the  expressions  is 

less  than  a factor  of  2 in  mass  at  a fixed  x . . 

cut 


It  should  be  noted  that  the  results  of  Figure  11  can  be 
displayed  in  a variety  of  ways,  each  answering  different  questions 
For  example,  the  amount  of  mass  ejected  up  to  time  t can  be 
obtained  by  replacing  xcut  with 


= a/RcZ  » a/  £(Z+1)  atj 


Z/(Z+1) 


and  Equation  (9)  can  be  verified.  Similarly,  the  amount  of  mass 
ejected  between  x1  and  x2,  or  R^  and  R2,  or  t^  and  t2  can  be 
obtained. 


Figure  11  Equation  (10)  results  compared  to  computed  JOHNIE  BOY  ejecta. 


The  experimental  measurements  of  the  JOHNIE  BOY  crater  are: 

R =21  meters  = apparent  crater  radius 

3i 

Rt  = 22  meters  = true  crater  radius 

V = 4900  m3  = ejecta  volume  (including  fallback  but 

e excluding  upthrust) 

The  ejecta  bulk  density  was  about  25  percent  less  than  that  of 

the  in  situ  soil.  Then  the  ejecta  volume  predicted  by  Equation  (10) 

is 


V 
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1.25 


ll  r3  (\  1 

3 rinal  L 2. 71+1 J 


0.5R 


3 

final 


= 4630  m3  for  R 

3L 

■ 5324  m3  for  RT 


where  is  in  good  agreement  with  the  measured  value  with  either 
choice  of  crater  radius. 


SECTION  3 


ENERGY  BALANCE  AND  FINAL  DIMENSIONS 


The  final  displacement  model  assumes  that  the  stream  line 
directions  do  not  change  when  the  flow  slows  and  that  energy 
balance  approximations  can  be  applied  independently  in  each 
stream  tube  to  stop  the  flow.  The  final  configuration  of  the 
mass  in  the  various  stream  tubes  then  defines  the  final  crater 
and  displacement  field.  This  process  is  shown  schematically 
in  Figure  12.  The  method  is  described  next. 

The  kinetic  energy,  KE,  the  gravitational  energy,  GE,  and 
the  total  distortional  energy,  DE,  in  the  steady  state  flow 
model  are  expressed  as  a function  of  Rc,  the  inner  radius  of 
the  mass  in  a stream  tube.  The  total  energy  in  the  stream 
tube  is  estimated  by  first  defining  the  condition 


KE(R+)  = GE(R*)  + DE(R*).  (11) 

c c c 

The  value  of  R+  is  obtained  from  Equation  (11)  and  is  used  to 
c 

estimate  the  total  energy,  EQ. 


= constant 


KE(R*)  + GE(R*)  + DE(R*). 


(12) 
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(c)  Stream  tube  2 flow  stops 
defining  another  final 
crater  *,ector 


(d)  Stream  tube  3 flow  stops, 
defining  crater  radius  and 
another  upthrust  point 


1 Igure  12  Schematic  configurations  of  mass  inside  various  stream 
tubes,  a,  b,  c,  and  d correspond  to  increasing  times. 


The  final  inner  radius,  Rf,  of  the  mass  in  the  tube  is  then 
evaluated  at  the  condition 


GE(Rf)  + DE(Rf)  = Eq.  (13) 

This  method  was  applied  to  two  stream  tubes:  the  axis  tube 

(in  which  case  Rf  is  the  final  crater  depth) , and  the  crater 
radius  tube  (in  which  case  Rf  = RAD) . The  latter  tube  was  de- 
fined by  the  condition  that  the  final  angle  0 equals  tt/2  con- 

c 

sistent  with  the  usual  definition  of  crater  radius. 

The  mathematical  development  associated  with  KE (R  ) , GE (R  ) , 

c c 

and  DE (R  ) is  given  in  Appendix  B and  only  selected  features  and 
c 

results  are  presented  here.  The  distortional  energy  assumptions 
are  the  most  severe,  and  require  a constant  yield  stress  approxi- 
mation to  obtain  a closed- form  solution. 

3.1  THE  CRATER  DEPTH  AND  RADIUS  MODEL 

The  Appendix  B model  predicts  the  crater  depth 


D = a1/Z  Y V(2« 

3L 


(14) 


where 
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The  corresponding  expressions  for  crater  radius  are 


RAD 


1/Z  1/ (2Z) 

a Yb 


(15) 


where 


Yb 


0.275 


YZ6 

(Z+l) 


r(i) 


0.2o  [h-(Z-2'/2]/(Z-1) 

2Z/(Z+1)  ' T.  >(2Z+l)/(Z+2) 


+ RAD 


TtSfr(r) 


and 


c = 0,  0.158,  0.268,  0.349,  and  0.41,  respectively 


for 


Z * 2,  2.5,  3,  3.5,  and  4 . 


Velocity  Cutoff — Cauchy  and  Froude  Scaling 

The  steady  state  velocities  that  existed  at  the  pc  ±tions  D 
and  RAD  before  strength  and  gravity  slow  the  flow  are 


D * a/DZ  « 1//7T 

A 


(16) 


RAD  * a/RADZ  - 1/ 


(17) 
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It  can  be  verified  that  y and  y.  have  no  dependency  on  a if 

cl  D 

the  gravity  term  can  be  neglected.  Then  Equations  (16)  and  (17) 
predict  that  all  cratering  flows  with  the  same  shape  (Z)  will 
freeze  at  the  same  cutoff  velocities.  This  is  consistent  with 
Cauchy  scaling. 

When  the  gravity  term  dominates,  then  the  cutoff  velocities 
are  proportional  to  the  square  root  of  the  depth  (or  radius) . 

This  is  consistent  with  Froude  scaling. 

3.2  CORRELATION  OF  DATA  BY  VELOCITY  CUTOFF 


Before  surface  burst  cratering  data  are  displayed  in  the 
velocity  cutoff  representation,  it  is  necessary  to  note  some 
qualifications  of  Table  1 cases. 


(a)  Case  1,  JOHNIE  BOY  has  reliable  experimental 
crater  data,  apparent  and  true.  The  computation 
employed  a weak  strength  model  with  a 1.6-bar  co- 
hesion that  vanished  in  the  regions  of  surface  spall. 
The  computed  crater  radius  was  satisfactory.  The 
computed  depth  was  about  40  percent  too  deep. 

(b)  Cases  2,  3,  and  4 had  reliable  experimental 
crater  data.  The  measured  value  of  the  plasticene 
yield  strength  at  atmospheric  pressure  was  about 
one  bar  (at  the  76°F  ambient  temperature  of  the 
experiments),  with  a large  error  due  to  creep  and 
a large  temperature  dependency.  The  calculations 
were  terminated  when  about  60  percent  of  the  energy 
coupling  to  the  ground  was  not  complete.  The  sub- 
sequent errors  in  a and  Z values  are  not  known. 

(c)  Case  5,  the  CACTUS  crater  dimensions  were 
measured  after  some  water  action  and  have  unknown 
reliability  except  that  the  depth  and  radius  measure- 
ments are  thought  to  correspond  respectively  to  minimum 
and  maximum  limits.  The  final  crater  was  not  computed. 


36 


(d)  Cases  7 and  8 were  very  similar  and  will  not  be 
distinguished.  There  is  no  experimental  data.  The 
computed  crater  dimensions  were  dominated  by  yield 
models  with  yield  strengths  that  were  very  large  com- 
pared to  the  successful  JOHNIE  BOY  model  in  the  low 
pressure,  highly  strained  region  of  the  cratering  flow. 


Cases  9 and  10  have  no  crater  data  and  Case  11  is  not  a 
surface  burst.  These  will  not  be  displayed. 


(e)  A dry  sand  cratering  experiment  (Z-Sl  of  Appendix  A, 
Table  Al)  was  conducted  in  the  scaled  geometry  of  Crse  3 
(plasticene  with  HOB  = 1) . The  flash  X-ray  photographs 
of  the  time-dependent  ejecta  lips  were  identical  in  the 
same  scale.  Therefore,  the  a,Z  values  of  Case  3 were 
scaled  to  the  size  of  the  sand  experiment  and  the  sand 
crater  depth  and  radius  are  included  as  data. 


The  cutoff  velocity  results  for  all  of  the  above  cases  are 
displayed  in  Figure  13.  The  various  lines  and  symbols  are 
explained  next. 


(a)  The  slanted  solid  lines  correspond  to  the  steady 
state  radial  velocities  using  the  volume  average  values  of 
ot  and  Z. 


R = a/RZ 

(b)  The  dashed  lines  labeled  D (Y)  and  RAD (Y)  are  the 
Equations  (16)  and  (17)  cutoff  velocities.as  a function 
of  Y in  bars  for  the  case  Z - 3#  p 55  2 . DiO)  and  RAD (0) 
are  *'he  gravity  controlled  limits. 
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R#  (neters) 


Figure  13  Velocity  cutoff  data  summary — see  text  for  symbols 


(C)  The  circle  and  square  data  points  correspond  to  apparent 
crater  depth  and  radius,  respectively.  The  labels  A,  T,  and 
C are  used  when  necessary  to  distinguish  the  apparent,  true, 
and  computer  craters. 

(d)  The  dotted  line  is  discussed  later. 


The  prediction  equations  were  inverted  to  define  the  effective 
yield  stresses  YD  and  Y^p  from  the  measured  crater  depth  and 
radius.  The  results  are  given  in  Table  3.  As  before,  the  symbols 
A,  T,  and  C are  used  when  necessary. 


TABLE  3 

EFFECTIVE  YIELD  STRESS  (BARS)  FOR  DEPTH  AND  RADIUS 


Case 

Y 

RAD 

CACTUS 

142.0 

-1.5 

ELK  76 

76.0  (C) 

87.0  (C) 

JOHNIE  BOY 

41.0  (A) 

1.6 (A) 

9.5  (T) 

0 . 3 (T) 

4.0  (C) 

0.6(C) 

HOB  = 0 

8.0 

18.0 

HOB  = 1 

3.4 

5.1 

HOB  = 2 

6.8 

4.5 

ZS-1 

4.1 

0.01 

Remarks 

Pacific  test  crater 
after  some  water  action 

Strong  model  of  layered  shales 

NTS  Area  18  alluvium 

Alluvium  model  with  1.6  bars 
cohesion 

Plasticene  clay  at  76°F 
Plasticene  clay  at  76°F 
Plasticene  clay  at  76°F 
Dry  sand 
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It  can  be  confirmed  that  the  6 plasticene  dimensions  can 
be  predicted  within  14  percent  by  a single  value  of  Y.  It 
appears  to  be  a nearly  ideal  plastic.  The  other  real  materials 
of  Table  3 are  not.  They  appear  to  have  very  little  strength 
with  respect  to  the  crater  radius  formation.  This  is  consistent 
with  the  moderately  successful  state-of-the-art  procedures  that 
predict  crater  radii  from  ballistic  trajectories  based  on  early 
computed  flows. 

3.3  ALTERNATIVE  GRAVITY  CONTROLLED  RADIUS  MODEL 

Before  describing  the  alternative  method,  it  is  useful  to 
note  the  features  of  the  usual  ballistic  method  which  is  applied 
to  the  flow  field  at  the  time  that  the  computer  crater  depth  is 
reached.  The  mass  points  are  assumed  to  be  non-interacting  as 
they  travel  along  ballistic  trajectories  into  mathematical  buckets 
which  are  fived  in  space,  usually  along  the  curved  surface  corres- 
ponding to  zero  vertical  velocity  at  crater  bottom  time.  The 
collected  mass  is  then  converted  into  volume  which  fills  the 
buckets  to  different  heights  to  define  the  final  crater  contour. 

The  moderate  success  of  this  method  prompted  the  development 
of  an  alternative  method,  intermediate  in  assumptions  between 
the  stream  tube  energy  balance  method  and  the  ballistic  method. 

In  this  new  method,  it  is  assumed  that  the  particles  proceed  on 
stream  ii*es  (rather  than  ballistic  paths)  but  that  spall  or 
some  other  effect  removes  the  particle  interactions  so  that 
energy  balance  applies  to  independent  particles  (rather  than 
the  entire  stream  tube).  The  final  crater  radius  is  then  de- 
fined by  the  particle  that  has  the  smallest  range  and  zero 
velocity  as  it  reaches  the  surface.  The  flow  model  provides 
a simple  way  of  selecting  this  particle. 
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The  locus  of  particles  with  purely  horizontal  velocities  lies 
along  the  line  cos  9^  = constant  = (Z-2)/Z-l)  in  the  first  order 
flow  model.  Consider  the  set  of  particles  i = 1,2, 3, .. .along 
this  line  with  coordinates  R.  ,0„  and  velocities  u..  The  kinetic 

th  l H p l 

energy  per  gram  at  the  i position  is  0.5  u^  . The  gravitational 

potential  is  gR.  cos  0„.  The  stream  line  intersects  the  surface 
^ “ Z— 2 

at  Rgi  given  by  (R^/R^)  = 1 - cos  0H«  It  can  be  shown  that 

only  one  of  the  particles  on  the  line  0„  has  the  property 

2 H 

gR.  cos  0„  = 0.5  u.  . The  crater  radius  prediction  is  then 
1 n 1 

h]1/(Z-2)  (18) 

with  the  constraint 

Rl  = u.2/(2g  cos  0H). 

Equation  (18)  can  be  expressed  in  a variety  of  convenient 
ways,  depending  on  the  form  of  the  velocity  field  data.  If  a ? i 
Z are  known,  then 


RAD  = R. 


cos  9 


RAD  = (Z-1)1/(Z'2)  [l+(Z-2)2/(2Z-3)]| 


1/ (2Z+1) 


(19) 


RAD  2q ( Z — 2 ) / ( Z — 1 ) / 1 \ 
^l+(Z-2)2/(2Z-3) j V2*1/ 


(2Z+1)/ ( Z— 2 ) 


(RAD) 


(20) 


The  velocity  cutoff  of  Equation  (20)  for  Z = 3 is  plotted  on 
Figure  13  as  the  dotted  line.  It  can  be  seen  that  this  method 
predicts  substantially  larger  crater  radii  than  the  previous 
method  ( RAD (0)  of  Figure  13].  For  ELK  76,  the  former  method 
predicts  a gravity  controlled  radius  of  about  115  meters, 
compared  to  about  175  meters  for  the  alternative  method. 
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SECTION  H 


ADVANCED  MODELS  AND  DATA 


The  cratering  flows  of  the  inspected  cases  could  be  approxi- 
mated with  fair  accuracy  by  the  first-order  flow  model.  This 
simple  model  provided  a framework  so  that  additional  assumptions 
could  be  included  to  predict  ejecta  properties  and  crater  dimen- 
sions. 

A variety  of  data  indicates  that  a second-order  model, 
allowing  for  spatial  variations  in  a and  Z,  would  reduce  the 
prediction  discrepancies.  Though  this  improved  order  model  has 
not  been  developed,  data  will  be  presented  to  indicate  reasonable 
features  that  should  be  included.  This  section  will  also  present 
miscellaneous  data  and  observations  which  are  not  understood  but 
which  may  be  significant. 

4.1  POROSITY 

It  was  thought  initially  that  porosity  effects  would  be 
severe  in  the  cratering  process.  Computed  and  experimental  data 
denies  this.  These  data  will  be  presented  next. 

Reference  2 reports  two  surface  burst  calculations  that  were 
identical  except  for  the  material  models.  The  first,  ELK  70 (L), 
employed  low  porosity  coral  layers  with  very  little  permanent 
compaction.  The  second,  ELK  75(H),  had  corresponding  layers 
with  high  porosity  and  substantial  permanent  compaction.  The 
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differences  in  the  two  shock  waves  were  large  but  the  stresses 
and  velocities  behind  the  shocks  were  very  similar.  The  horizontal 
velocities  versus  time  and  stresses  versus  time  at  a common  point 
are  displayed  in  Figures  14  and  15.  The  large  differences  in  the 
shock  waves  (due  to  porosity)  are  not  observed  in  the  steady  state 
flows  that  begin  at  about  5 msec. 

The  displacement  fields  in  the  cratering  regions  were  also 
similar  at  12  msec  when  the  calculations  were  terminated.  This 
can  be  seen  by  comparing  the  grid  distortions  displayed  in  Fig- 
ures 16  and  17.  The  current  crater  depth  and  radius  of  the  low 
porosity  case  are  larger  than  those  of  the  high  porosity  case 
by  1 percent  and  3 percent,  respectively. 

These  computation  observations  are  supported  by  flash  X-ray 
photographs  of  the  time-dependent  crater  lip  growth  seen  in  the 
plasticene  and  sand  experiments  that  had  identical  sources 
(tangent  spheres  of  HE) . The  results  are  displayed  as  the  circle 
and  square  points  in  Figure  A-23  of  Appendix  A.  The  differences 
are  small  in  spite  of  the  large  differences  in  porosity. 

Apparently  the  cratering  flows  are  controlled  primarily  by 
the  source  parameters,  with  material  properties  becoming  important 
only  at  late  times  when  strength  effects  appear. 

The  fact  that  the  crater  lip  growths  are  the  same  for  the 
non-porous  plasticene  and  the  very  porous  dry  sand  also  has 
implications  with  respect  to  the  conjecture  that  the  airblast, 
seeping  into  the  porous  material,  will  change  the  ejection  rate. 
Experimental  high  explosive  airblast  data,  scaled  to  the  size 
corresponding  to  Figure  A-23  enters  the  negative  pressure  phase 
at  45  psec  at  the  ranges  of  the  ejecta.  If  air  is  entrained  in 
the  sand  at  high  pressure,  it  does  not  appear  to  affect  the 
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ejection  rate  compared  to  that  of  the  corresponding  plasticene 
within  the  scatter  of  the  data. 

4.2  PLASTICENE  DISPLACEMENTS  AND  VARIABLE  Z 

Several  of  the  cratering  experiments  of  Appendix  A had 
color  banding  in  the  medium  to  establish  the  final  displacement 
fields.  These  can  be  seen  in  Figures  A-5  through  A-8  of  Appen- 
dix A.  The  points  shown  on  Figure  A-5  for  a plasticene  case 
were  compared  to  the  Equations  (3)  and  (4)  expressions 

RZ+1  . Z+l  = Rz+1  (3) 

o c 

(R/R  )Z_2  = (1-cos  6) /(1-cos  6J  (4) 

o o 


with  Z and  Rc  being  unknown.  Trial  and  error  fits  were  performed 

and  it  was  found  that  the  best  fit  values  of  Z and  R_  were 

c 

essentially  constant  at  given  value  of  6 when  Equation  (3)  was 
employed.  The  value  of  R was  found  to  agree  well  with  the 

C 

measured  crater  radius  at  this  angle. 

The  fits  for  the  14  points  on  the  axis  are  displayed  in 
Figure  18a.  It  can  be  seen  that  the  best  fits  are 

Z(0°>  = 1.9  t 0.1 

R (0°)  = 11.2  ± 0.1  cm 

c 

The  measured  crater  depth  was  11.1  t 0.1  cm  which  is  in  good 
agreement. 


U8 


5 


Orig.  slant  range.  RQ  (ml 


Figure  18  Fits  to  Equation  (3).  (a)  Top,  14  axis  points 

(b)  Bottom,  11  surface  pairs 


The  22  surface  points  were  grouped  in  11  pairs  to  eliminate 
asymmetries.  The  fits  are  shown  in  Figure  18b.  The  results  are 

Z (90°)  = 3.0  + 0.2 

R (90°)  = 12.25  ± 0.1  cm 

c 

The  value  of  Rc  is  in  good  agreement  with  a crater  radius  con- 
structed by  a smooth  extrapolation  through  the  lip  spall  region. 

The  6 points  at  60°  were  fit  to  Equation  (3)  as  individual 
points  (not  in  pairs) . The  results  were: 

Z(9  = 60°)  = 2.9  i 0.1 

Rc(6  = 60°)  = 12.0  ± 0.1  cm 

Again  the  fit  value  of  R agreed  well  with  the  experimental  radius 

c 

(slant  range)  of  the  crater  wall  at  this  angle. 

Equation  (4)  was  also  applied  to  these  60°  points  and  the 
results  were  surprising.  Here  the  particles  were  grouped  into 
pairs  to  show  the  tendencies  more  clearly. 

Inner  pair:  Z(60°)  = 1.8  t 0.1,  Rc(60°)  = 10.8  t 0.1  cm 

Middle  pair:  Z(60°)  = 1.9  ± 0.1,  Rc(60°)  = 9.8  i 0.1  cm 

Outer  pair:  Z(60°)  = 2.0  t 0.1,  R (60°)  = 9.4  i 0.1  cm 

c 

These  last  results  are  not  understood.  The  other  cases  (Fig- 
ure A- 6 through  A- 8)  have  not  been  analyzed. 
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4.3  COOPER  DISPLACEMENTS 


The  final  radial  displacements  of  the  previous  example  were 
consistent  with  the  first  order  radial  displacement  model.  This 
suggests  a connection  between  the  model  and  Cooper* s displacement 
correlations  that  were  observed  in  both  soil  and  rock  (see 
Reference  1 for  rock  example).  Cooper* s correlations  are: 


d„ 

s 

K„  V4/3 

Y -3 

H 

H 

o 

dv 

K V4/3 

Y "3 

V 

o 

for  2 

< 

Yq/v1/3 

S 6 

where  d„  and  d are  respectively  the  horizontal  (outward)  and 
H v 

vertical  (upward)  displacements  at  the  horizontal  range  Y . 

V is  the  apparent  crater  volume.  and  Kv  were  found  to  vary 
with  gauge  depth  and  explosive  DOB.  The  values  KH  = Ky  = 0.45 
gave  reasonable  agreement  with  surface  gauge  data  for  above  the 
surface  bursts. 

The  first  order  model  can  be  used  to  construct  expressions 
for  dH  and  dy.  Ignoring  higher  order  terms  in  d/YQ,  the  result 
is 

du  = R 2+1  Y ”2  (sin  6)Z+1  11  + (Z-2)  cos  6/(1  + cos  9)|/(Z+1) 

n C O 

dv  = RCZ+1  Y0"Z  (sin  Q)Z  I (Z-2)  (1-cos  0)  - cos  0|/(Z+1) 

where  R is  the  inner  radius  of  the  displaced  mass  in  a stream 
c 

tube  that  passes  through  the  gauge  at  Yq  and  tf.  The  model  ex- 
pressions have  Cooper's  form  if  Z = 3 and  R is  proportional  to 
vl/3. 


51 


However,  Cooper's  correlations  are  for  sets  of  gauges  at  the 

same  depth  rather  than  the  same  angle  and  a direct  comparison  of 

the  methods  is  possible  only  for  near  surface  gauges  where 

0 = tt/2.  In  this  case  the  model  results  for  Z = 3 is 

4 -3 

d„  = d = 0.25  RAD  Y where  RAD  is  the  displacement  crater 

H . V ° 1/3 

radius  at  the  surface.  A relationship  between  RAD  and  V is 

required  in  order  to  proceed. 


Cooper  (Reference  1)  reported  that  the  MINERAL  ROCK  and 

MINE  ORE  apparent  crater  radii  R were  related  to  the  apparent 

1/3  A 

crater  volume  by  R.  - 1,2  V ' . If  R.  is  equated  to  the  displace- 

A A 4/3  -3 

ment  radius  RAD,  then  the  model  result  is  d„  = d = 0.52  V ' Y 

H v o 

which  is  16  percent  larger  than  Cooper's  result.  Exact  agree- 
ment is  obtained  if  RAD  = 1.16 


It  is  concluded  that  Cooper  scaling  and  the  radial  dis- 
placement model  are  compatible.  It  is  suggested  that  the 
displacement  crater  can  be  approximated  by  the  apparent  crater, 
rather  than  the  true  crater. 


4.4  Z VALUES  AND  VERTICAL  MOMENTUM 


Computed  and  experimental  data  show  that  the  cratering  and 
ejection  process  occurs  on  a very  long  time  scale  compared  to 
that  of  the  direct  and  airblast  loading  of  the  cratering  region. 
This  implies  that  the  steady  state  flow  field  is  a coasting  pro- 
cess with  zero  total  vertical  momentum  (up  plus  down) . This 
places  a constraint  on  allowable  values  of  Z.  At  early  times 
before  surface  uplift  occurs,  it  is  easy  to  show  that  the  only 
constant  value  of  Z that  gives  zero  vertical  momentum  is  Z = 3. 

This  value  is  also  the  only  solution  of  Equations  (1)  and 
(2)  that  gives  irrotational  flow.  The  significance  of  this  is 
not  known  at  this  time. 


Vertical  momentum  data  were  obtained  in  the  flow  fields  of 
the  10  surface  burst  cases  of  Table  1 as  part  of  the  output  of 
the  ZEE  code  (described  in  Section  4.5).  This  code  processed 
printout  data  grouped  as  concentric  surfaces  about  ground  zero 
(see  the  dashed  lines  of  Figure  lc  for  an  example) . The 
vertical  momenta  of  the  points  on  a given  surface  were  summed, 
starting  from  the  axis  point,  to  obtain  the  accumulated  vertical 
momentum  as  a function  of  0 along  each  surface.  More  than  100 
surfaces  were  processed  in  this  manner,  selected  from  the  10 
surface  burst  cases  of  Table  1.  Several  different  times  were 
used  in  two  of  the  cases.  Though  this  large  amount  of  data  has 
not  been  thoroughly  analyzed,  it  was  seen  that  the  total  vertical 
momentum  contained  between  0 and  90°  was  nearly  zero  in  all  cases. 
Apparently,  the  vertical  momentum  is  nearly  zero  in  concentric 
shells,  as  well  as  over  the  entire  flow.  This  implies  that  the 
angular  dependency  of  Z is  such  that 
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[/ 


Z(0)R  sin  0 d0/ 


TT  / 2 

/ 


R sin  0 d0 


fixed  R 


4.5  THE  ZEE  CODE  AND  FLOW  CHARACTERIZATION 

A code  called  ZEE  was  written  and  used  to  process  printout 
data  to  examine  flow  parameters  in  detail.  The  ZEE  inputs  are 
point  values  of  coordinates,  velocities,  and  densities.  The 
outputs  are  values  of  candidate  expressions  and  graphical 
displays. 

The  code  first  constructs  point  values  of  Z from  Equation  (2) 
and  then  finds  the  corresponding  a from  Equation  (1) 
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2 + R0(l+cos  0 ) / (R  sin  0) 


a 


R R 


Z 


The  Z values  are  noisy  at  small  0 because  of  the  sin  0 denominator. 
This  should  be  allowed  for  in  the  contour  plot  shown  later. 

The  units  of  a depend  on  the  value  of  Z,  and  a will  be  noisy 
if  Z varies  even  slightly.  It  is  more  meaningful  to  measure  flow 
strength  at  a point  in  terms  of  an  equivalent  crater  dimension 
or  velocity  cutoff  in  a standard  material.  This  was  done  in 
several  ways.  One  method  combined  Equations  (14)  and  (15)  in  a 
crude  way  to  define 


R1  = [“2  yJ1/(2z>  jcos2e  + [l+  ( Z— 2 ) 2]  sin2  ej 

with  Y = 1 bar  and  g = 0.  A contour  plot  of  R^  will  be  shown 
later. 


Numerous  other  expressions  were  processed.  The  vertical 
momentum  and  its  properties  were  described  earner . Another 
useful  expression  that  will  be  displayed  is 


ZBAR (R, 0 ) 


0 

/ 


Z (R, 0)  sin  0 d0/  / sin 


0 

/ 


0 d0 

■'fixed  R 
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This  quantity,  evaluated  at  0 = 90°  and  averaged  over  R should  be 
similar  to  the  volume  average  Z of  the  first-order  model. 


Another  promising  quantity  is 

sin  0 d0 I 

fixed  R 

This  quantity,  evaluated  at  0 = 90° , should  be  similar  to  the 
volume  average  R used  in  Figure  3. 

The  Z,  ZBAR,  R^,  and  R representations  described  above  are 
displayed  in  Figures  19  through  22  for  the  same  flow  field.  The 
case  is  ELK  66  (tangent  HE  sphere)  at  30  ysec.  The  shock  wave  is 
just  beyond  the  range  of  the  data.  The  shock  tail  extends  inward 
to  about  4.5  cm  and  accounts  for  the  upward  curvatures  in  this 
region. 


This  set  of  four  displays  was  obtained  for  all  of  the  10 
surface-burst  cases  of  Table  1.  Several  times  were  used  in  two 
of  the  cases  to  confirm  that  the  time  dependencies  were  small. 
Although  it  is  not  obvious  in  the  early  time  data  of  Figures  19 
through  22,  the  other  data  showed  that  the  variations  with  R were 
small  in  the  region  well  behind  the  shock.  The  0 dependency  re- 
mained, however. 


Except  for  spatial  scale,  it  was  seen  that  the  patterns  of 
the  10  cases  were  similar  and  could  be  grouped  in  order  of  simil- 
arity. The  half-buried  HE  sphere  (ELK  67)  was  at  one  extreme, 
and  the  HE  sphere  elevated  by  one  radius  (ELK  69)  was  at  the  other. 
The  other  patterns  were  intermediate  with  the  tangent  HE  sphere 
(ELK  66)  being  most  similar  to  the  5-Mt  surface  burst  (ELK  76). 
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The  data  has  not  been  analyzed  in  sufficient  detail  to 
uncover  a quantitative  means  of  describing  and  ranking  these 
four.  It  is  felt  that  this  can  be  done  and  that  a refined  flow 
model  would  result. 
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SECTION  5 


SUMMARY  AND  CONCLUSION 


5 . 1 SUMMARY 

It  was  found  that  the  cratering  flows  observed  in  the  ten 
computed  surface  burst  cases*  could  be  approximated  by  a first- 
order  model  to  predict  time-dependent  trajectories,  displacements, 
and  ejecta.  An  elastic-perfectly  plastic  material  assumption 
together  with  an  energy  balance  approximation  was  used  to  predict 
final  crater  sizes  and  shapes  and  to  explore  the  scaling  rules 
and  the  roles  played  by  material  strength  and  gravity. 

Data  for  extending  and  improving  the  models  were  collected 
and  processed  but  not  adequately  analyzed.  Some  general  features 
of  an  improved  description  were  established  though  the  final  model 
was  not. 

5.2  CONCLUSION 

The  approximate  agreement  of  the  existing  models  to  a variety 
of  cases  and  applications  is  the  most  signifcant  result  of  this 
study.  It  is  concluded  that  a reasonable  framework  has  been 
established  that  unifies  the  processes  of  cratering,  close-in 
displacements  and  ejecta,  and  provides  a basis  for  improving 
understanding. 

The  MINE  ORE  computed  flow  was  also  examined.  The  ground  shock 
from  the  high  explosive  source  was  not  strong  enough  to  initiate 
a significant  region  of  steady  state  flow  due  to  an  unrealistically 
strong  granite  model. 
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APPENDIX  A 

EXPERIMENTAL  REPORT  OF 

CRATERING,  DISPLACEMENT,  AND  EJECTA  PROCESSES 
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Various  experimental  tests  were  performed  using  specially 
prepared  test  beds  and  small  charges  to  examine  the  dynamics 
of  crater  development  and  the  permanent  post-shot  deformations 
that  occur  as  a result  of  crater  formation  under  controlled 
laboratory  conditions. 

Two  materials  that  simulate  earth  media  were  selected  for 
the  experimental  test  bed.  Plasticene  clay,  a suspension  of 
white  marble  dust  in  a petrolatum  and  wax  matrix  with  a density 
of  1.69  gm/cm3,  was  selected  to  simulate  a clay  earth  media. 
Nevada  Type  47  sand  was  selected  to  simulate  sandy,  alluvial 
soil.  This  sand  is  a fine-grain,  high-purity,  quartz  sand  with 
an  uncompacted  density  of  1.53  gm/cm3.  For  experiments  per- 
formed in  this  program,  the  sand  was  kept  dry  and  in  an  uncom- 
pacted state.  Separate  test  beds  were  made  for  the  two  test 
materials.  The  test  beds  were  pits  sunk  in  the  ground  with  the 
top  surface  level  with  the  surrounding  soil.  These  pits  were 
8-feet  square  by  4-feet  deep  and  were  enclosed  in  a 2-inch-deep 
concrete  apron  at  the  surface  level. 

As  tests  were  performed  in  the  plasticene  material,  it 
became  evident  that  crater  size  was  strongly  dependent  on  clay 
temperature.  Attempts  were  made  to  characterize  the  dependency 
of  plasticene  strength  on  temperature  variations.  Unconfined 
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compressive  strength,  shear  strength,  and  ball  drop  penetration 
tests  were  performed  in  an  attempt  to  obtain  a clearer  view  of 
plasticene  behavior.  Unfortunately,  the  data  were  difficult  to 
interpret  because  of  creep.  Plasticene  strength  is  apparently 
extremely  strain-rate  and  temperature  dependent  and  cannot  be 
well  characterized  in  an  essentially  static  test,  such  as  an 
unconfined  compressive  strength  test  or  a shear  strength  test. 
The  ball  drop  data,  being  less  affected  by  creep,  appears  more 
promising. 

1.1  PERMANENT  CLOSE-IN  DEFORMATION  AND  DISPLACEMENT  TESTS 

The  first  experimental  test  series  examined  the  permanent, 
close-in  deformations  and  displacements  that  occur  as  a result 
of  crater  formation.  Surface  displacements  were  measured  by 
placing  wooden  pegs  (1/16-inch  in  diameter  by  approximately 
3/4-inch  long)  as  markers  in  the  clay  or  sand  and  measuring 
their  pre-  and  post-shot  distance  from  a permanent  bench  mark. 
The  centerline  of  the  shot,  relative  to  the  bench  mark,  was  also 
measured  before  the  shot,  and  thus  the  radial  displacement  could 
be  obtained.  Crater  profile  and  volume  were  obtained  by  using  a 
crater  contour  gauge.  Since  the  craters  were  symmetrical,  a 
profile  taken  through  the  centerline  of  the  crater  was  suffi- 
cient to  determine  crater  volume.  Figure  A1  is  a photograph  of 
the  crater  contour  gauge  over  a symmetrical  sand  crater.  Since 
the  sand  craters  were  relatively  large,  the  contour  gauge  was 
limited  to  taking  a profile  along  a radius  from  the  center. 

Close-in  crater  displacement  and  deformation  were  obtained 
by  placing  a series  of  multicolored  columns  vertically  beneath 
the  centerline  of  the  shot.  In  the  plasticene  media,  the 
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columns  were  prefabricated  into  a slab  of  plasticene  (Figure  A2)  , 
which  was  placed  into  the  test  bed.  Molten  clay  was  then  poured 
around  the  slab  to  ensure  a voidless  fit  (Figure  A3) , and  the 
explosive  was  placed  on  the  centerline  of  the  slab.  Following 
the  shot,  the  slab  containing  the  columns  was  removed  from  the 
test  bed  and  was  sectioned  along  the  centerline  to  obtain  a 
pictorial  representation  of  the  close-in,  permanent  deformation 
(see  Figures  A4  through  A8) . 

The  close-in  deformation  of  the  sand  test  bed  was  obtained 
by  placing  multicolored  columns  of  epoxy- loaded  sand  along  the 
centerline  of  the  shot.  The  sand  surrounding  the  columns  (and 
the  columns  themselves)  was  loaded  with  1 percent,  thermosetting, 
powdered  epoxy.  After  the  test,  an  aluminum  box  was  forced 
gently  around  the  area  containing  the  conditioned  sand  and  the 
crater  centerline.  This  box  and  the  sand  in  it  were  then  heated 
until  the  epoxy  was  cured,  and  then  the  competent  sand  was 
removed  for  sectioning  and  inspection  (see  Figures  A9  through  A12) . 

To  determine  whether  or  not  the  1-percent  epoxy  loading 
alters  the  crater  formation  processes,  a shot  was  fired  in  the 
sand  test  bed,  but  with  no  epoxy  in  the  sand.  The  results  from 
this  test  showed  that  a loading  of  1 percent,  powdered  epoxy  in 
the  sand  has  no  effect  on  crater  formation.  A comparison  of  the 
crater  profiles  is  shown  in  Figure  A13)  . 

Surface  displacement  records  were  taken  by  placing  wooden 
pegs  (approximately  1/16-inch  diameter  by  3/4-inch  length)  as 
markers  and  measuring  their  pre-  and  post-shot  distance  relative 
to  a permanent  bench  mark.  These  pegs  were  protected  from  move- 
ment caused  by  the  airblast  effects  by  a large  steel  and  wood 
cover  that  was  placed  over  the  row  of  pegs,  almost  but  not  quite 
touching  the  sand.  Crater  profile  and  volume  were  obtained  by 
using  a crater  contour  gauge. 
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£u  - 


Figure  A2 


Multicolored  columns  of  plasticene 
made  into  a slab  for  placement  in 


before  being 
the  test  bed. 


Figure  A3 


Columns  in  place  in  test  bed  with  2 -inch-diameter 
sphere  of  C-4  ready  for  firing. 


Figure  A- 5 Experimental  data,  shot  Z-C2,  2-in.-diam  C-4  sphere, 
tangent  to  surface  of  clay.  Clay  temperature  = 76  F. 

The  circled  points  are  analyzed  in  Section  4 of  main  text 
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Figure  A-12  Sand  crater,  shot  Z-52  (tangent  sphere). 


Crater  radius,  cm  fromq 


Figure  A-13  Effects  of  epoxy  loading  on  crater  formation  (pure 

sand  versus  sand  loaded  with  1-percent  powdered  epoxy) 
Crater  formed  with  a 2-inch-diameter  sphere  of  C-4 
tangent  to  surface  of  sand. 
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Surface,  permanent  displacements  are  shown  plotted  against 
crater  radius  in  Figures  A14  through  A18.  For  these  plots,  both 
displacement  and  original  radius  were  divided  by  crater  volume 
to  the  one- third  power.  The  straight  line  through  the  data  is  a 
fit  to  the  results  obtained  by  Cooper. 

It  can  be  seen  that  the  half-buried,  spherical,  explosive 
cases  of  Figures  A14  and  A15  (corresponding  to  shots  Z-C5  and  Z-C4 

at  temperatures  of  63.5°F  and  90°F)  give  results  that  are  in 
fair  agreement  with  Cooper's  data.  The  surface  tangent  cases  of 
Figures  A16  and  A18  also  yield  fair  agreement.  The  elevated  shot- 

case  of  Figure  A17  shows  a much  different  range  dependency  than 
that  of  the  surface  burst  data  of  Cooper,  a result  that  is  not 
surprising. 

The  plasticene  crater  experiments  corresponding  to  scaled- 
up  versions  of  the  calculations  ELK  67,  66,  and  69  employed 
high-explosive  charge  diameters  of  1.375,  2.0,  and  2.5  inches, 
respectively,  to  obtain  crater  sizes  that  utilized  the  test  bed 
efficiently.  The  corresponding  scale  factors  of  length  to 
relate  the  ELK  experimental  charge  sizes  to  the  respective  ELK 
calculations  are  0.573,  0.394,  and  0.315  before  detonator  cor- 
rection. When  the  type  RP-2  detonator,  its  volume,  and  its 
high-explosive  content  are  taken  into  account,  the  previous 
scale  factors  are  changed  slightly  to  become  0.579,  0.3945,  and 
0.3153,  respectively.  The  results,  including  the  results  of  the 
sand  experiment  which  were  not  calculated,  are  shown  in  Table  Al. 

Figure  A19  displays  two  experimental  crater  sections  for 
each  of  the  scaled  experimental  data  of  ZC-6,  zc-5,  and  ZC-7  for 
the  76°F  cases  corresponding  to  the  calculations  ELK  67,  66,  and 
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Figure 


A-14  Shot  Z-C5,  surface  horizontal  displacement  versus 

range  (Loth  scaled  by  cube  root  of  crater  volume;  in 
plasticere  1-3/8  in.  diam  C-4  sphere,  half  buried. 


1/3 


0.01 


0.001 


Shot  Z-C2,  surface  horizontal  displacement  versus 
range  (both  scaled  by  cube  root  of  crater  volume)  in 
plasticene  2-in.-diam  surface  tangent,  C-4  sphere. 


1/3 


Figure  A-17  Shot  Z-C3,  surface  horizontal  displacement  versus 
range  (both  scaled  by  cube  root  of  crater  volume) 
plasticene  elevated  1 radius. 
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(d)  Includes  small  correction  for  detonator 

(e)  Scaled  crater  displaced  mass  M * = V ' x a (where  c = 1.69  gm/cm3  for  clay  and  1.53  gm/cm3  for  sand) 
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69.  The  dashed  circles  show  the  nearly  spherical  shapes.  The 
differences  in  the  crater  sections  of  the  same  experiment  suggest 
an  experimental  nonreproducibility  of  about  0.1  cm  in  the 
"scaled"  data. 

1.2  CRATER  EJECTA  GROWTH  TESTS 

A second  test  series  was  performed  to  examine  the  growth  of 
the  crater  ejecta  "splash"  and  the  corresponding  crater  lip. 
Independently  triggered,  30  kV  flash  X-ray  units  were  used  to 
obtain  time-resolved  radiographs  of  the  growth  of  the  crater  lip. 
Since  only  three  X-ray  units  were  available,  the  time-resolved 
series  is  necessarily  constructed  from  successive  tests.  Care 
was  taken  to  insure  that  the  test  bed  was  in  the  same  configura- 
tion for  each  successive  test.  Fresh  sand  and  clay  were  used  for 
each  shot,  and  the  temperature  of  the  sand  was  monitored  to 
insure  the  equivalent  temperature  for  each  shot. 

Figures  A20,  A21,  and  A22  show  the  sequential,  time-resolved 
radiographs  of  the  growth  of  the  ejecta.  Two  test  series  were 
conducted  in  the  plasticene  test  bed;  one  series  was  conducted  in 
the  sand  test  bed.  The  wires  observable  in  the  radiographs  are 
electrical  leads  to  the  detonators. 

Figure  A23  is  a graph  of  the  minimum  diameter  of  the  ejecta 
"splash"  versus  time.  The  values  plotted  have  been  scaled  down 
to  match  the  ELK  66,  67,  and  69  calculations.  Table  A2  tabulates 
the  results  of  this  experimental  series,  as  well  as  the  scaled- 
down  values.  It  is  interesting  to  note  that  the  curves  of  the 
ejecta  "splash"  for  the  two  tangent  tests  fall  quite  closely  to 
each  other,  ever,  though  the  material  properties  are  significantly 
different.  This  observation  implies  a scalability  of  crater 
growth  for  scaled  HOB  geometry  before  strength  effects  appear  at 
later  times.  This  point  is  supported  by  computational  data  that 
will  be  reported  later. 
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2-inch  diameter  sphere  of  C-4  explosive 
plasticene  clay 


IT 


Figure  Acl  Sequential  flash  radiographs  of  crater  ejecta 
formation  in  plasticine. 
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1-inch  diameter  sphere  of  C-4  explosive,  tai 
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25  usee  after  ZT  50  usee  after  ZT 


1.3  CRATER  EJECTA  MASS  TESTS 


The  final  experimental  test  series  was  conducted  to  examine 
the  actual  mass  of  ejecta  excavated  from  the  test  media  during 
explosive  cratering.  To  accomplish  this,  a method  was  developed 
that  allowed  an  18-inch-diameter  plasticene  hemisphere  of  known 
weight  to  be  lowered  into  a prepared  receptacle  in  the  plasticene 
test  bed.  The  shot  was  fired,  and  the  hemisphere  with  crater  was 
then  removed  and  reweighed  to  determine  mass  loss.  Figure  A24 
shows  the  receptacle  ready  for  the  hemisphere  of  plasticene. 

Figure  A25  shows  the  plug  in  place  and  a half-buried  shot  ready 
for  firing.  The  results  of  this  test  series  are  shown  in  Table  A3. 


TABLE  A3 

RESULTS  OF  EJECTA  MASS  EXPERIMENTS 


Shot 

Clay 

Temperature 

(°p) 

Configuration 

Crater 

Volume 

(cm3) 

Eject 

Mass 

(gm) 

CRADY-19 

60 

2-in.  C-4  sphere 
tangent  to 
surface 

1509 

951 

CRADY-20 

55 

1-3/8-in.  C-4 
sphere  1/2 
buried 

1164 

1121 

C RADY- 21 

56 

2-1/2-in.  C-4 
sphere  elevated 
1-1/4-in.  above 
surface  of  clay 

815 

57 

Three  plasticene  cratering  experiments  to  measure  ejected 
mass  were  relevant  to  the  data  base  for  prediction  purposes.  The 
results  scaled  down  to  the  ELK  66,  67,  and  69  calculations  size 
are  shown  in  Table  A4. 
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Figure  A24  Hemispherical  receptacle  in  plasticene  test  bed 
coated  with  a thin  layer  of  Adiprene  L-100  as  a 
mold  release. 


Figure  A25  Hemispherical  clay  in  place  in  receptacle,  1-3/8-inch- 
diameter  C-4  sphere  half  buried  in  plasticene  and  ready 
for  shot. 


TABLE  A4 


SCALED 

PLASTICENE 

EJECTA  RESULTS 

Scaled  Quantity 

Shot 

Crady-20 

Shot 

Crady-19 

Shot 

Crady- 

HOB  (cm) 

0 

1 

2 

Crater  radius  (cm) 

5.3 

4.0 

2.4 

Crater  depth  (cm) 

4.7 

3.1 

1.1 

Crater  volume  (cm3) 

226 

92.6 

25.5 

Ejecta  mass  (grams) 

218 

58.4 

1.79 

Lip  mass  (grams) 

121 

43.5 

8.90 

Clay  temperature  (°F) 

55 

60 

56 

:t  is  unfortunate  that  the  three  shots  of  Table  A4  were  not 
at  the  same  temperature.  However,  the  combined  data  of  Tables  A1 
and  A4  afford  a means  of  performing  reasonable  interpolations  to 
relate  data  at  a common  temperature  (hence  at  a common  material 
strength)  for  comparisons  to  the  computed  data  base. 
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APPENDIX  B 

STREAM  TUBE  ENERGY  BALANCE  EQUATIONS 
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MATERIAL  STRENGTH  MODEL 


Although  the  strain  field  is  prescribed  by  the  steady  state 
flow  model,  the  distortional  energies  also  depend  on  the  material 
model.  It  was  necessary  to  approximate  real  materials  with  a 
constant  shear  modulus,  y,  and  a constant  von  Mises  yield  stress, 
Y,  in  order  to  proceed. 

Typical  values  of  Y/y  vary  from  about  0.012  to  0.0013  going 
from  aluminum,  through  Westerly  granite  and  mild  steel,  to 
plasticene  clay.  Therefore,  the  approximation  Y/y  < 1 was  used 
to  eliminate  small  terms  in  various  expressions. 

The  yield  stress,  Y,  was  used  to  limit  the  second  invariant 
of  stress  deviation  /JjJ  in  the  plastic  domain.  Although  the 
stress  deviators  were  different  in  the  axis  and  radius  stream 
tubes,  the  values  of  /3jJ  differed  only  a few  percent  at  corres- 
ponding points.  Therefore,  the  axis  tube  expressions  for  the 
plastic-elastic  transition  radius,  R* , and  the  plastic  and  elas- 
tic distortional  energy  densities,  np  and  ri£,  can  be  applied  in 
both  tubes  with  little  error. 


PLASTIC-ELASTIC  TRANSITION  RADIUS,  R* 

The  plastic-elastic  transition  radius,  R* , was  expressed  as 


c 

[Rc/R*]Z+1  = 1 - eXP 

(z  + M 

l 3Z  / 

Applying  Y/y  < 1,  the  result  is 


L 3*  \ ii 

\Z + 1/  Y 


1/(Z+1) 


(B-l) 


which  has  values  ranging  from  about  6 to  12  for  Z = 2 (contained 
bursts)  and  from  about  4 to  6 for  Z = 3 (typical  surface  burst 
average  value) . 


DISTORT I ON AL  ENERGY  DENSITIES  AND  INTEGRALS 


The  plastic  and  elastic  distortional  energy  densities  np 
and  nE  can  be  evaluated  when  Y and  y are  constant.  Using  the 
approximation  Y/y  < 1,  the  results  are 


nF(R  < R*)  = YV(6y) 


np(R  < R*)  = 


np  (R*  < R)  = 


- y2/(3)j)  - imr  1 - (VR)Z+1 


(m)  7 (VR) 


2(Z+1) 


np(R*  < R)  = 0 


The  stream  tube  volume  element  is  A R dR  where  A_  is  the 

o o 

cross  sectional  area  at  unit  range.  The  distortional  energies 

contained  between  R,  and  R0  in  a stream  tube  with  A = 1 are 

12  o 


J VR  < R*)RZdR  = — 

Rx  E 2 (Z  + 1) 


Z + 1\  Y nZ+l  _ DZ+1 

“Tz  ) U R2  ■ Rl 


/ 1,, 


< R)R  dR  = 


2 (Z  + 1) 


3Z  \ u 
TT/  Y 


,2 (Z+l) 


np(R  < R,)R  dR  = 


5 k1  - «rl 


YZ  „Z+1 
? Rc 

(Z  + ir 
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AXIS  TUBE  KINETIC,  GRAVITATIONAL,  AND  DISTORTIONAL  ENERGIES 

The  kinetic,  gravitational,  and  distortional  energies  in  the 
axis  stream  tube  are  given  by  integrals  of  the  form 

00 

A Qf  f (R)  RZdR 
Rc 

The  results  are  given  below. 


KEaxis/Ao 

= R1_Z 

pa2 

(B-2) 

Rc 

2(Z  - 1) 

GEaxis/Ao 

- RZ  + 2 

■ Rc 

pg  x, 

(B-3) 

(Z  + 2) 

( i 

DEaxis/Ao 

■-r1 

YZ  * ’i  + ln  f 3Z  v]( 

(B-4) 

TT7T71 1 

* RZ+1 

YZ  6 

(B-4a) 

Kc 

(z  + ir 

The  logarithm  term  of  Eq.  (B-4)  varies  by  about  +20  percent  for 
common  materials.  It  is  replaced  by  the  average  value  of  5 to 
obtain  the  Eq.  (B-4a)  approximation. 

RADIUS  TUBE  KINETIC  ENERGY 

The  radius  tube  kinetic  energy  is  the  same  as  that  of  the 
axis  tube  with  correction  terms  involving  Rc/RAD  that  cannot  be 
evaluated  unless  Eqs.  (11),  (12),  and  (13)  of  Section  3,  main 
text  are  iterated.  However,  since  KE^p  is  used  only  in  Eq.  (11), 
it  is  sufficient  to  express  the  correction  terms  as  functions  of 
rVrad.  It  can  be  shown  that 
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^,1/fZ+l) 

w 


< rVrad  < 

w 


1/ (Z+2) 


with  the  lower  limit  corresponding  to  GE  < DE  and  the  converse 
applying  at  the  upper  limit. 


The  correction  terms  were  evaluated  at  these  limits  for 
values  of  Z = 2,  3,  and  4 and  it  was  found  that  KE^p  (with 
R - R+)  could  be  approximated  by 

KErad  ~ °-4  KEaxis  t1  + <Z  - 2)2]  (B‘5) 

for  application  in  Eq.  (11)  with  a few  percent  error. 


RADIUS  TUBE  GRAVITATIONAL  ENERGY 


The  change  in  gravitational  potential  in  the  radius  tube 

also  contains  terms  involving  R /RAD  with  values  R+/RAD  applying 

c c 

in  Eq.  (11)  and  R /RAD  = 1 applying  in  Eq.  (13).  Paralleling  the 

w 

KErad  Procedure'  the  results  are 


where 


gerad 

GEaxis  f RAD 

for 

Eq. 

(13) 

gerad 

= GE  f 

axis  R+ 

for 

Eq. 

(11) 

(B-6) 

(B-7) 
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Values  of  are  0,  0.186,  0.315,  0.410,  and  0.482,  respectively 

for  Z = 2,  2.5,  3,  3.5,  and  4.  The  error  in  Eq.  (B-8)  is  a few 
percent,  rising  to  about  10  percent  in  the  Eq.  (B-9)  approximation. 

RADIUS  TUBE  DISTORTIONAL  ENERGY 

The  plastic-elastic  transition  radius  is  well  beyond  the 

outer  radius  of  the  mass  when  R+  < in  the  radius  tube.  As 

c c 

before,  the  results  are  the  same  as  those  of  the  axis  tube  with 
a correction  that  depends  on  the  value  of  Rc/RAD.  Using  previous 
procedures,  the  results  are 


derad  “ °-23 

DE 

axis 

for 

Eq. 

(13) 

(B-10) 

DERAD  = °-32 

DE 

axis 

for 

Eq. 

(ID 

(B-ll) 

with  a few  percent  error. 

INTERIM  SUMMARY 

Let  the  subscripts  i = 1,2  apply  respectively  to  the  axis 
and  radius  stream  tubes  and  define  the  quantities 
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n 
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Equations  (11) , (12) , and  (13)  give 


Ki 


-t + — r 

Di  + Rc  Gi 


1/(2Z) 


E_  = 2K. 


Z-l 


|1/(Z+1) 


R, 


Di  + Rf  Gi 


(B-12) 


(B-13) 


(B-14 ) 


where  iterations  are  required  to  evaluate  Eqs.  (B-12)  and  (B-14) 
if  the  gravity  term  and  strength  term  are  comparable. 

Equation  (B-14)  can  be  evaluated  directly  by  employing  *-he 

i i 

approximations  <Pi>  = 0.5(Di  + D^)  and  <G^>  = 0.5^  + G^)  and 
writing  Eq.  (B-14)  as 


77  77 

(Rc/Rf)  <Pi>+(Rc/Rf)  Rf-Gi> 


1/(Z+1) 


The  appropriate  estimates  of  R*/Rj  are  (jj 
respectively  for  the  <D^>  and  <G^>  terms. 


1/(Z+1) 


and 


j) 


1/  (Z  + 2) 


The  result  is 


Ri 


Ki 


, /TTTzT(zTY) 7.  .('Sz+U/JzVJJ 

L (?)  <Di>  + r£  (j)  -v 


1 1/ ( 2Z ) 

(B-l 5) 
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